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' Abstract 

• 

The semirelativistic Hamiltonian H = (3\Jm 2 + p 2 + V(r), where V(r) is a central 
potential in 9ft 3 , is concave in p 2 and convex in p = \pp 2 ~. This fact enables 
us to obtain complementary energy bounds for the discrete spectrum of H. By 
extending the notion of 'kinetic potential' we are able to find general energy bounds 



on the ground-state energy E corresponding to potentials with the form V = 
Yli a if^( r )- I* 1 the case of sums of powers and the log potential, where V(r) = 
J2 q =£o a(Q)sgn(q)r q + a(0) ln(r), the bounds can all be expressed in the semi-classical 
form 

E « min I pJm 2 + \ + a(q)sgn(q)(rP(q)) q + a(0) ln(rP(0)) > . 
{ r J 

'Upper' and 'lower' P-numbers are provided for q = —1,1,2, and for the log 

potential q = 0. Some specific examples are discussed, to show the quality of the 

bounds. 



PACS: 03.65. Ge, 03.65. Pm, ll.10.St 
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1. Introduction 

The Hamiltonian H = K + V for the problems we study has the feature that 
either the kinetic energy K or the potential energy V is nonlocal. The most im- 
portant example is the semirelativistic "spinless-Salpeter" Hamiltonian [1-4] given 
by 

H = K(p) + V(r) = v 7 ™ 2 + P 2 + V(r), (1.1) 

where r = ||r||, r G 9? 3 , and p = ||p|| = y/p*. In this form at least K(p) is non- 
local in configuration space and is defined as a multiplicative operator in momentum 
space. That is to say, Kip is defined to be what we get when tjj is transformed to 
momentum space, the multiplicative operator K is applied, and the result is trans- 
formed back to coordinate space. Nonlocality is the main source of difficulty for this 
class of problems. We study the discrete spectra of these Hamiltonians by the use of 
approaches that make use of convexity and of spectral information already obtained 
concerning related problems. In an earlier paper [4] we studied the relationship be- 
tween Ho = K + h and H = K + goh, where g(h(r)) is a smooth transformation 
of a 'base' potential h(r). For cases in which g had definite convexity, one could 
then employ the so-called 'kinetic-potential' [5] formalism and 'envelope theory' [6] 
to construct upper or lower bounds to the discrete eigenvalues of H by using the 
known spectrum of Hq. 

The present paper has two distinct aspects: we turn our attention firstly to 
the convexity of K, as a function of p or p 2 ; and then we look at potentials 
that are a sum of terms V = J2 i In Section 2 we extend the kinetic-potential 
formalism to include more general kinetic-energy operators than the Schrodinger 
form K = p 2 studied earlier [5]. We have already found [7] some implications of 
the fact that K is concave in p 2 . The spectral implications of the convexity of 
K in p demand a new analysis. By using Jensen's inequality [8] we are able to 
construct a framework in Section 3 which accommodates both cases. Although it 
may not be immediately apparent, the treatment of potential sums also leads to 
an interesting convexity analysis of a completely different type [9, 10]. This in turn 
yields an optimized general lower bound for the bottom of the spectrum expressed 
in terms of the kinetic potentials generated by 'component' problems K + V^K We 
discuss this in general terms in Section 4. In the special case in which the component 



Convexity and potential sums for Salpeter-like Hamiltonians 



page 3 



problems are pure-power potentials or the log potential, that is to say 

V(r) = a(q)sgn(q)r q + a(0) ln(r), (1.2) 

the lower bound (Section 5) and also upper bounds obtained by variational methods 
(Section 6) can all be expressed in terms of a semi-classical expression of the form 



E ~ min < J m 2 + + ^ a(q)sga(q)(rP(q)) q + a(0) ln(rP(0)) 



(1.3) 



It is the goal of this paper to develop a general theory which leads to such a result 
and, in particular, to determine P-numbers which guarantee that the approximation 
(1.3) is an upper or lower bound. We shall find the appropriate upper and lower 
P(q) for the cases q = —1, 1,2, and for the log potential q = 0. In Section 7 we 
apply our general results to some specific examples. 

Although we obtain very concrete results in the end, our study begins with a 
somewhat abstract viewpoint. We now make a few general remarks that will help 
motivate these starting considerations. An idea that runs through the work is one 
well known to those who study non-linear problems: we try to use transformations 
to make the most of any soluble problem that is at hand, or, at least, one for which 
we have a good approximation. The setting for our ideas is geometrical. We sup- 
pose that we have an exact solution (or good bounds) for a 'base' problem with 
Hamiltonian ap 2 + (3h(r), and we are interested in a Hamiltonian of the form 
H = k{p 2 ) + g(h(r)), where k and g are monotone increasing smooth transfor- 
mations. It follows that the 'tangent spaces' to H are Schrodinger operators with 
the general form = a + bp 2 + ch(r), where the parameters {a, b, c} depend 

on the contact vector t. What we look for is a theory that would allow us to 
deduce spectral information about H from the 'known' spectrum of its tangents 
_H"(*) . For example, if the transformation functions k and g are both concave, we 
would expect to obtain upper bounds via the spectral inequality H < Mu- 
tatis mutandis, a complimentary theory is possible with operator tangent spaces of 
the form = a + bp + ch(r), where K{p) is convex in p and g(h) is convex 

in h : this leads to energy lower bounds. Since we have already explored [4] the 
potential transformation g, the main thrust of the present paper concerns kinetic- 
energy transformations with 'base problems' respectively p 2 + h and p + h and to 
applications of the results when the potential in H is a sum of terms. The more 
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general case in which neither transformation k(p) nor g(h) is the identity follows 
immediately by combining the present and earlier results; hence we do not need 
to discuss this natural generalization in detail here. Coulomb 'components' present 
a special difficulty for a lower bound with kinetic energy p because the operator 
p — v/r has no discrete eigenvalues. However, for our main concern, the Salpeter 
Hamiltonian, we are able to make use of the very good Coulomb lower bound of 
Martin and Roy [11] and so incorporate the Coulomb contribution smoothly into 
our general formulation. 

2. Variety of kinetic potentials 

The discrete eigenvalues of the self-adjoint operators we study may be char- 
acterized variationally. Thus the bottom of the spectrum E of H = K + V is 
given by E = inf (tfj, Htfj), where the infimum is taken over all normalized functions 
in the domain T>(H) C L 2 (K 3 ). The idea behind kinetic potentials is to perform 
the minimization in two stages: we first find the constrained minimum V (K; s) of 
(ip, Vip), keeping the mean kinetic energy (ifj^Kif;) = s constant; then, we recover 
E by minimizing over the kinetic energy s > 0. Thus we have 

V(K;s)= inf (ip,Vip) E = minis + V(K;s)\ . (2.1) 

4>eT>(H) s>o 1 

(V,V)=i 

We call the function V(K; s) the kinetic potential of V associated with the kinetic- 
energy operator K; we shall write simply V(s), if the kinetic-energy operator K 
is fixed or is clear from the context. It follows immediately from the definition 
that the kinetic potentials absorb a positive coupling parameter in the sense that 
cV(s) = cV(s). We note also that the elementary 
Comparison Theorem 

V (1 \s)<V (2 \s) =► E^<EW 

follows immediately from (2.1). The arguments we use are not restricted to dimen- 
sion N = 3 : this choice allows us to illustrate the general results with some explicit 
well-known examples, without the distraction of the operator dependencies on N. 

The reason for using this description of the spectral problem is that it lends it- 
self to some interesting approximations. Firstly, we have shown in the Schrodinger [5] 
and Salpeter [4] cases that, if V(r) = g(h(r)) and g is monotone increasing and has 
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definite convexity, then the approximation V(s) ~ g(h(s)) leads to lower bounds if 
g is convex and upper bounds if g is concave. In Section 4 of the present paper we 
shall extend to general K the result obtained earlier [10] for the Schrodinger case 
that kinetic potentials are subadditive, that is to say 

V(r) = h {1) (r) + h {2 \r) V(s) > h (1) (s) + h {2) (s). (2.2) 

The lower energy bound then immediately follows from the above-mentioned com- 
parison theorem for kinetic potentials. There is more to this result than meets the 
eye: it generates the optimum of a family of lower bounds; the details will be given 
in Section 4 below. Extensions to sums with more than two terms (or, further, 
to mixtures generated by an integral) are immediate. The principal limitation is 
that each potential term alone, when added to the kinetic energy, must, for large 
enough coupling, support a discrete eigenvalue. Thus V(r) = —1/r + r is allowed 
but V(r) = 1/r + r is not. 

The 'component' kinetic potential h(s) = h(K; s) can be constructed by use of 
a Legendre transformation from the eigenvalue function E = F(v), in which F(v) 
is the bottom of the spectrum of H = K + vh(r), as a function of the coupling v. 
In the Schrodinger case H = p 2 + vh(r) we have shown that F(v) is concave [5] 
and moreover the kinetic potential for K = p 2 is given in terms of F(v) by the 
transformation 

8 = F(v) -vF'(v), h(K;s) = F'(v). (2.3) 

The concavity of the eigenvalue function F(v) has been proved for the Schrodinger 
case [5] and the Salpeter case [4] by the application of a simple variational argu- 
ment. By exactly similar reasoning we can show that the eigenvalue function for 
the operator p + vh(r) is also concave in v. Moreover, the Legendre transfor- 
mation (2.3) is generic: it is valid for all kinetic-energy operators K. This is an 
immediate consequence of the concavity of F(v), as the following equations clearly 
demonstrate: 

F(v) = mm{F(u) - uF'(u) + vF'(u)} = mm{s + vh(K; s)}. 

u>0 s>0 

Our principal assumption concerning K is that it is at once a convex function of p 
and a concave function of p 2 . This convexity is clearly true for our most important 
example, the relativistic kinetic energy K = ^/m 2 + p 2 ; however, we shall use this 
specific form only when we need to. 
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We now turn from the general to some very specific results. We shall need to 
have at our disposal some 'component' kinetic potentials for the operators K + 
v sgn(g)r g , where K = p , or the Schrodinger case K = p 2 . By elementary scaling 
arguments we can show that the dependence of the energy functions on the coupling 
v are given by 



and 



p + v sgn(q)r q =>- E = F^(q; v) = F^\q; l)v^ (2.4a) 
p 2 + v sgn(g)r 9 E = F (2) (q; v) = F (2) (q; l)v& . (2.46) 



The Legendre transformation F <-> h given above in (2.3) now allows us to deduce 
the precise forms of the corresponding kinetic potentials. For convenience we choose 
to write the kinetic potentials so obtained in a special way. We change variables for 
the mean kinetic energy s in the two cases respectively to s = 1/r and s = 1/r 2 . 
It then follows from (2.3) by straightforward algebraic computations that the kinetic 
potentials for h(r) = sgn(q)r q have similar convenient forms, namely 



h(p; 1/r) = sgn(q)(P {1) (q)r) q 



(2.5a) 



and 



h(p 2 ; 1/r 2 ) = sgn(q)(P (2 \q)r) q , (2.56) 

where the P-numbers are defined in terms of the v = 1 eigenvalues E^ l \q) = 
F( l \q; 1), i = 1, 2, respectively by the explicit formulas 



K = p 



E^(q) 



l + q 



q\ , Q > -!> Q o 



(2.6a) 



and 



K = p z 



£( 2 )(g) 


2 + q 


q 


l + q/2 




2 



q> -2, q^O. 



(2.66) 



The energies are related to the kinetic potentials by specific realizations of the 
general formula (2.1): for example, we have in this spectral representation 

p + v sgn(q)r q =3- E = mm | - + v sgn(g) (^P^(q)r^ |, q > — 1, q ^ 0. 

J (2-7) 
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One of our side goals is purely esthetic, namely we wish to end up with 'attractive 
formulas': after the changes of variable from s to r, the kinetic potentials look like 
the original power potentials themselves, but with the P-factors inserted. We turn 
now to the base potential h(r) = ln(r) and find by scaling arguments that 

p + vln(r) => F (1 \v) =vF {1 \l) -vln(v) (2.8a) 

and 

p 2 + v\n(r) F (2 \v) = vF {2 \l) - jUln(v). (2.86) 
Consequently we obtain from the transformation (2.3) 

h(p; 1/r) = ln(pW(0)r), h(p 2 ; 1/r 2 ) = ln(p( 2 )(0)r), (2.9) 

where 

P (1) (0) =exp(pW(0)-l), P^ 2) (0) = -J=exp ^ 2 )(0) - . (2.10) 

For the discussion of examples we shall need to have some specific P values. For 
the cases q = — 1, 0, 1, 2 we supply some of these numerical values in Table 1. This 
table has an eigenvalue symmetry because of the operator equivalence p+r 2 ~ p 2 +r; 
it also has two omissions corresponding to q = — 1, because p — 1/r has no discrete 
eigenvalues. We offer now a solution to this Coulomb difficulty. As we shall make 
clear in Section 4, viable Coulomb P-numbers are needed for lower bounds. For our 
most important application K = \J p 2 + m 2 , a lower bound to the bottom of the 
spectrum of H = K — v/r is provided by the Martin-Roy bound [11] 

*> ei (») = m(i±4^\ v <\. (2.11) 

The condition v < \ is a little more restrictive than the fundamental operator 
restriction v < 2/ir : it was proved by Herbst [12] that a Friedrichs extension exists 
for H only if the Coulomb coupling is sufficiently small. The Coulomb lower bound 
has the same scaling law with respect to m as does the exact energy: although 
m originates in the Hamiltonian inside the square root of the kinetic-energy term, 
it appears in the eigenvalue and in its lower approximation simply as an overall 
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factor [4]. Now we construct a independent P-representation for this lower bound. 
We write (as a definition of Pl(v) ) 



e L (v) = min { Jm 2 + - — — - . (2.12) 




An elementary calculation then shows that (2.11) and (2.12) imply 
Pl(v) = eL,(v)/m. This serendipitous discovery fills the gaps in Table 1, 
and will allow us to include the Coulomb component in our lower-bound energy 
formula for sums of potential terms: we must make the substitution 

v v mv 1 , n _ „. 

v < -. (2.13) 



p(i)(_l) r P L (y)r e L (vy 2 



3. Complementary convexity: p + V and p 2 + V 

The principal result of this section is best expressed in terms of kinetic poten- 
tials by the following 

Theorem 1. If E is the bottom of the spectrum of the Hamiltonian H = K + V, 
and the kinetic- energy operator K is at once convex in p and concave in p 2 , then 
it follows that 

min {K(s) + V(p;s)} < E < mm {K(s) + V(p 2 ; s 2 )} . (3.1) 

It makes sense here to speak of K(p) as though p were a real variable since, by 
definition, the action of the operator K is effected via the Fourier transform. We 
shall now prove this result by an application of Jensen's inequality [8] and kinetic 
potentials defined in (2.1). We consider first the left-hand inequality of the theorem. 
If if) is a normalized function in the domain T>(H) of H = K + V, then, since 
K = K(p) is convex in p, by Jensen's inequality, we have 

E= inf AU>,K(p)if>) + U>,V1>)}> inf |K((V>,W0) + V^)} . 
(V,V)=i (V-,V-)=i 

That is to say, 

E > min inf {K((ib, pip)) + (ib, Vib)} = min \K(s) + V(p: s)\ . 
~ s>0 4,eT>(H) J s>0 1 Vjr /J 



Convexity and potential sums for Salpeter-like Hamiltonians 



page 9 



The proof of the upper-bound inequality is very similar: we write K(p) = k(p 2 ), 
where k is concave; then, setting t = s 2 , we arrive at the inequality 

E < min {k(t) + V(p 2 ;t)} =mm{K(s) + V(p 2 ;s 2 )}, 

which establishes the theorem. □ 
This result is an essential ingredient in the proof of the sum approximation in the 
next section. We now look at an example, namely the Salpeter problem with a linear 
potential. We have 



H = ^/m 2 +p 2 + V(r), V(r) = vh(r) = vr, (3.2) 

where v is a positive coupling parameter. In terms of the convenient variable r > 
the two kinetic potentials from (2.5) are 

h{p; 1/r) = pW(l)r, h(p 2 ; 1/r 2 ) = p( 2 )(l)r, (3.3) 

where the P- numbers are provided in Table 1. Theorem 1 then immediately yields 
the bounds 



mm 

r>0 



|v / m 2 + r- 2 + vP (1) (l)r} < E < mm | ^m 2 + r~ 2 + vP {2) (l)r} . (3.4) 

In Figure 1 we plot these bounds as a function of m for the case v = 1 . If we 
combine Theorem 1 here with Theorem 2 of Ref. [4] (to the effect that g o h > g o h 
when g is convex) we obtain the following class of examples. We suppose that V(r) 
is monotone increasing and convex in h(r) = r then the two theorems together yield 
the lower bound 

E > mm j^m 2 + + wy(P (1) (l)r)| . (3.5) 

Of course, if V(r) is concave, then we get an upper bound by the same expression 
provided we use p( 2 )(l). It is perhaps important to note that with P = P^ 2 \r) 
an upper bound would be obtained for every choice of r in the expression on the 
right-hand side; the expression in (3.5) is however only a lower bound a priori at 
the minimum point. 

4. The sum approximation: lower bounds 

Since further generalization easily follows, we first look at the problem of the 
sum of two potential terms. We assume that each potential vh^ (r) alone, when 
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added to the kinetic-energy operator K, has a discrete eigenvalue at the bottom of 
the spectrum for sufficiently large 'coupling' v. We express our result in terms of 
kinetic potentials and prove the following 

Theorem 2 If E is the bottom of the spectrum of the Hamiltonian H = K + V, 
and the potential V is the sum V(r) = h^(r) + h^ 2 \r), then it follows that the 
sum of the component kinetic potentials yields a lower bound to V, that is to say 

V(K;s)>h (1) (K;s) + h (2 \K;s). (4.1) 

We shall now prove this theorem, which is in effect an optimized Weyl lower 
bound [13-15]; this remark will be clarified below, after the proof of the theorem. 
From the definition (2.1) of kinetic potentials we have 

V(K;s)= inf (^,V^) = mf U, + h {2) ) V>V 
(V,V)=i (i^)=i 

(tP,K4>)=s (ip,Kip)=s 

But the latter minimum mean-value is clearly bounded below by the sum of the 
separate minima. Thus we have 

V(K; s) > inf U, h (1) ifj) + inf U, h (2) t(j) = h (1) (K; s) + h (2) (K; s) , 

(V,V)=i W»,V0=i 
(4>,K4>)=s (4>,K4>)=s 

which inequality establishes the theorem. □ 

Another approach, which would eventually yield an alternative proof of the 
theorem, exhibits the relationship between Theorem 2 and the classical Weyl lower 
bound [13-15] for the eigenvalues of the sum of two operators. Let us suppose 
that ^ is the exact normalized lowest eigenfunction of H = K + V", so that 
H^f = E^>. If the positive real parameter w is bounded by 1, < w < 1, then 
E = (K + V)^) may be written as follows: 

E = w (^K + i- /i (1) (r)^ ^ + (1 - w) (ty, + h {2 \r)^j * 

> w inf (ij), ( K + — /i (1) (r) ) ip) 
(V,^)=i 

+ (1-w) inf (V, (k+ — ?— /i (2) (r)W 



Convexity and potential sums for Salpeter-like Hamiltonians 



page 11 



That is to say, in terms of component kinetic potentials, we arrive at Weyl's in- 
equality for the lowest eigenvalue of the sum H = wK + + (1 - w)K + : 

E > w min { s H (K: s) \ + (1 — w) min < s H (K: s) 

~ s>0 { W J s>0 I 1-W 

Since w is an essentially free parameter in the last expression, we may optimize the 
Weyl lower bound with respect to the choice of w and this forces the individual 
values of s at the minima, {si(w), s 2 (w)}, to be related. More specifically we find 
from the individual minimizations over s, 

E > £(w) = W8l (w) + (1 - w)s 2 (w) + h {1) (K; 8l (w)) +h {2) {K-s 2 {w)), 



where 



dh (1) dh {2) 
w = -— (K;si(w)), and 1 - w = -— (K;s 2 (w)). 

OS OS 



The critical condition £'{w) = for the subsequent maximization of the lower 
bound over w then yields s\(w) = S2(w). Thus the best lower energy bound is 
given by 

E > mm | s + h {1) (K; s) + h (2) (K; s) } . 

The kinetic-potential inequality of Theorem 2 leads, of course, to the same energy 
lower bound: the optimization just performed above is therefore seen to be auto- 
matically 'built in' by the formalism. 

It follows immediately from the above kinetic-potential comparison theorem 
and coupling-parameter absorption that a lower bound to the lowest energy E of 
the Hamiltonian H = K + £\ Cih^\r), {q > 0}, is provided by the formula 

E > mm js + ^c^ W (if;s) j . (4.2) 

Similarly we can extend this result to 'continuous sums' such as V(r) = 

This general theory becomes practically useful when we have good information 
concerning the components. More particularly, we must have some exact component 
kinetic potentials, or lower bounds to them. Outside the well-explored Schrodinger 
case K = p 2 , such analytical results are rather sparse. We look at the interesting 
class of power-law potentials in the next section. 
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5. Sums of powers and the log potential 

For power-law potentials and the relativistic kinetic energy K = a/ m 2 + p 2 we 
have discussed some lower bounds in Section 3 and we shall now turn these to our 
advantage. The link between Theorem 1 and Theorem 2 derives from the observation 
that the equation K(s) = \Jm 2 + s 2 allows us to change the minimization variable 
s — > r = 1/s. In the first stage of minimization, we have used Jensen's inequality for 
the lower bound (see proof of Theorem 1); this allows us to keep (ip,pip) = s = 1/r 
constant at first, and then later minimize over s, or, equivalently, over r. We can 
also easily accommodate a further positive kinetic-energy parameter (3. Thus we 
immediately arrive at 

Theorem 3 A lower bound to the lowest eigenvalue of the semirelativistic spinless- 
Salpeter operator 

H = f3\J m 2 + p 2 + a(q)sgn(q)r q + a(0) ln(r), 
where (3 > and the potential coefficients a(q) > are not all zero, is given by 



E > minJ/3Wm 2 + 4 + 5^ a (?) s g n (?)( p(1) (9) r ) 9 + a(0)ln(P (1) (0)r) >, (5.1) 
r> ° { r q ^o J 

where, for the Coulomb component q = — 1, we make the substitution 

<*(-!) . 0*™ 0v ( 2 y q(-l) 1 

P(!)(-l)r e L (v)r r \l + JY^W ) ' P 2' 1 ' J 

The problem presented for the lower bound by the fact that p — v/r has no discrete 
spectrum was discussed in Section 2. We have no simple P-number P(— 1) but we 
could derive a 'running' P (2.13) from the Martin-Roy energy bound (2.11); the 
positive factor (3 has been inserted in (5.2) by elementary scaling. We shall look 
at applications of Theorem 3 in Section 7 when we also have at our disposal the 
upper-bound P-numbers derived in Section 6. 

6. Variational upper bounds 

The lower bound for sums discussed in the previous two sections has the attrac- 
tive feature that if the component kinetic potentials are exact and only one term is 
present, then the result is exact. We are unable to construct a general upper bound 
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with this feature. Instead we use a trial wave function <fi = cexp(— \ar v ) with a 
scale parameter a > and two other parameters {c, v}, and we apply this wave 
function to the entire problem. One degree of freedom c is used to guarantee nor- 
malization, and the scale parameter a > is expressed in terms of a new variable 
t > chosen in such a way that the scale minimization is of an expression with the 
same form as the lower bound. Initially we use here t rather than r since, during 
the discussion, we shall need to refer to the potential function V(r). The choice of 
the remaining parameter v > is left for later optimization. 

If we suppose that c has already been chosen so that \\<j)\\ = 1, and, for 
computational convenience, we use Jensen's inequality, we then obtain the following 
upper energy bound: 



E <S = /V™ 2 + (<P, P 2 <P)+ y>, |^a(o)sgn(Q)r 9 + a(O)ln(r)|0j . (6.1) 
Now, for each fixed v > 0, we define a new scale variable t > by the following: 
(^) = (0, -A0) = a>/» {^y^^sf 1 - \- (6-2) 

Using this definition of t , we can go on to define the 'upper' P-numbers V(u, q) 
by the relations 

i r(2±3) 

far"*) = aq/v ^ =(V(u,q)t) q , o^O, (6.3a) 

and 

(0,ln(r)0) = \n(V(v,0)t). (6.36) 

If we now rename the scale variable t = r, and minimize the upper bound £ with 
respect to scale, we arrive at 

Theorem 4 For each v > 0, an upper bound to the lowest eigenvalue E of the 
Salpeter operator 

H = 13^/m 2 +p 2 + ^ a(q)sgn(q)r q + a(0) ln(r), 
where (3 > and the potential coefficients a(g) > are not all zero, is given by 
E < min < pJm 2 + \ + a(g)sgn(a)(P(z/, Q)r) 9 + a(0) ln(P(z/, 0)r) I , (6.4) 
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where the upper F '-numbers are provided by the formulas 

•—WW- «" 

arid ^ is the dig amma function ip(t) = F'(t)/T(t). 

Apart from the special Coulomb considerations pertaining to the lower bound (5.1), 
that formula is essentially identical to the upper bound (6.4): we simply have to use 
the correct P-numbers in each case. 

7. Examples 

We have now assembled the P-numbers for our energy-bound formulas (5.1) 
and (6.4). We shall use the lower P-numbers in Table 1, lower 'running' P- formula 
for the Coulomb component (2.13), and the formulas (6.5) for the V(u,q) corre- 
sponding to the variational upper bound (6.4). The class of problems we are thus 
immediately able to consider have the following explicit Hamiltonian form: 

H = I3\fm 2 +p 2 - a/r + bln(r) +cr + dr 2 , a,b,c,d>0, (7.1) 

where (3 > 0, and the potential parameters {a, 6, c, d} are not all zero. We look 
at two examples. In the first, illustrated in Figure 2, we look at the linear-plus- 
Coulomb potential V(r) = — 0.1/r + 0.25r and compare the energy bounds {L, U} 
we find, as functions of the mass m, with some very accurate numerical values 
(center curve) obtained by minimizing the expectation value of the Hamiltonian 
in a 25-dimensional trial space. In the next graph, Figure 3, we plot the energy 
bounds alone, for the same potential and a wider range of values of the mass m. 
As a second example we consider the broad linear combination V(r) = — 0.1 /r + 
0.25 ln(r) + 0.25r + 0.25r 2 and plot in Figure 4 the energy bounds as functions of 
the mass. These illustrations give a clear indication of the quality of the bounds 
that the theory yields. 

8. Conclusion 

The principal theoretical results of this paper are the complementary bounds of 
Theorem 1, and the sum-approximation lower bound, Theorem 2. In order to arrive 
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at these results we needed first to extend the notion of 'kinetic potential' to allow 
for more general kinetic-energy operators than the Schrodinger form K = p 2 . The 
complementary bounds are based on the assumption that K is a convex function 
of p and also a concave function of p 2 , assumptions clearly satisfied by our prime 
example and principal motivation, the relativistic kinetic energy K = j3^J m 2 + p 2 . 
The inequality of Jensen then allows us to learn approximately how special mean 
values of the problem, the eigenvalues of H, depend on the operator parameters. 

By combining Theorem 1 of this paper with Theorem 2 of our earlier paper [4] 
we obtain a general theory applicable to 'operator manifolds' of the form H = 
K(p) + g(h) with, on the one hand, tangent spaces spanned by the Schrodinger 
operators ap 2 + bh(r) + c, and, on the other, by complementary operators of the 
form ap + bh(r) + c. Given the correct convexities of K and g, energy bounds 
immediately follow. We looked at one example of this type of problem near the end 
of Section 3; and the results were exhibited in Figure 1. 

A completely different lower bound is provided by Theorem 2, which may be 
thought of as a spectral expression of the sum structure of the potential, namely 
the subadditivity of the corresponding kinetic potential, as a sum of components. 

In order to make practical use of these theoretical results we need some definite 
spectral information about component problems. This is provided by the family 
of pure-power potentials V(r) = sgn(q)r q . For this family we are able to take 
advantage of known eigenvalues, or bounds to them, and of simple upper bounds 
obtained with the aid of Jensen's inequality and a two-parameter family of trial 
functions. All our component results can then be expressed in terms of certain 
P-numbers (or, for the lower Coulomb case, q = — 1 , by a P- function), which 
are required by the general lower- and upper-bound formulas of Theorems 3 and 
4. These formulas illustrate the effectiveness of the theoretical results and provide 
recipes for approximate solutions to an interesting class of semirelativistic spectral 
problems. 
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Table 1. Eigenvalues for v = 1 and corresponding P- numbers [given by (2.6) 
and (2.10)] for the Coulomb, log, linear and harmonic-oscillator potentials. The 
eigenvalues have been computed numerically and are rounded so that the E^\q) 
are lower bounds and the E^ 2 \q) are upper bounds (and similarly for the derived 
P-numbers). The Coulomb lower bound is treated differently because H = p — 1/r 
has no discrete spectrum. 
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Complementary upper ( U ) and lower ( L ) bounds (3.4) on the lowest eigenvalue 
E(m) of H = a/ m 2 + p 2 + r plotted against m . 
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Figure 2. 

Lower bounds ( L ) by (5.1) and upper bounds ( U ) by (6.4) for the lowest eigenvalue 
E{m) of H = \/ m 2 + p 2 — 0.1/r + 0.25r plotted against m . The upper bound ( U ) 
used the wave- function parameter v = 1.6. The central curve is a very accurate 
upper bound found by a variational exploration in a 25-dimensional trial space. 
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Figure 3. 

Lower bounds ( L ) by (5.1) and upper bounds ( U ) by (6.4) for the lowest eigenvalue 
E(m) of H = \/ m 2 + p 2 — 0.1/r + 0.25r plotted against m : this is a continuation 
of the graph in Figure 2 to larger m. 
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Figure 4. 

Lower bounds ( L ) by (5.1) and upper bounds ( U ) by (6.4) for the lowest eigenvalue 
E(m) of H = a/tb 2 + p 2 - 0.1/r + 0.25 ln(r) + 0.25r + 0.25r 2 plotted against m. 
The upper bound ( Z7 ) used the wave- function parameter z/ = 1.4. 



